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Abstract 

The  author,  in  his  masters  thesis  [3],  constructed  a  catalog  of  442 
Bernoulli  sampling  plans  which  approximately  minimize  the  maximum  expected 
sample  size  among  all  plans  which  guarantee  certain  O.C.  probability  require¬ 
ments.  Fifty-two  of  these  plans  (which  would  appear  to  be  of  greatest  prac¬ 
tical  interest)  are  presented  in  this  report.  A.S.N.  curve  comparisons  are 
made  with  plans  based  on  the  Wald  sequential  probability  ratio  test  and  the 
fixed  sample  size  test  which  guarantee  the  same  O.C.  probability  requirements. 
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1.  Introduction 


In  a  paper  by  D.  Freeman  and  L.  Weiss  [1],  Weiss  proposed  using,  for  a 
Bernoulli  population,  a  sequential  sampling  plan  with  an  A.S.N.  curve  that 
has  the  smallest  maximum  among  all  sampling  plans  satisfying 

P„(plan  accepts  the  population)  >  1-a  if  t>  <, 0  , 

(1.1) 

Pe(plan  accepts  the  population)  8  if  0  :>6j 

where  6  represents  the  percent  defective  in  a  population  and  O<0_j<6j<1,  a>0, 
6>0,  a*B<l.  He  presented  the  theory  for  constructing  plans  which  approxi¬ 
mately  minimize  the  maximum  expected  sample  size,  and  with  Freeman  developed 
a  few  examples  of  such  plans.  For  all  of  these  Weiss  plans,  the  number  of 
observations  required  by  a  plan  is  a  bounded  chance  variable;  the  fact  that 
the  bound  is  known  for  each  plan  is  of  obvious  practical  importance.  In  Sec¬ 
tion  2  the  theory  behind  the  Weiss  test  is  discussed. 

Other  important  work  in  the  area  of  sequential  sampling  has  been  done 
by  A.  Wald  [9]  and  the  Statistical  Research  Group  at  Columbia  University 
(S.R.G.)  [7,8].  Wald  introduced  sequential  sampling  plans  which  minimize  the 
A.S.N.  curve  at  0_j  and  0^  among  all  plans  satisfying  (1.1),  and  the  S.R.G. 
developed  and  cataloged  sequential  plans  which  in  many  instances  (but  not 
always)  possess  A.S.N.  curves  that  are  lower  than  those  of  single-  or  double¬ 
sampling  plans  which  satisfy  (1.1). 

In  a  masters  thesis  [3],  the  author  presented  a  catalog  of  442  sequen¬ 
tial-sampling  plans  as  suggested  by  Weiss.  The  present  paper  will  display  52 
of  these  plans  which  are  appropriate  for  small  values  of  0_j,0j  and  thus  are 
considered  to  be  of  greatest  practical  interest.  The  A.S.N.  curve  character¬ 
istics  of  these  plans  will  be  compared  wj*h  those  of  corresponding  Wald  and 
fixed  sample  size  plans. 
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2.  Theoretical  Development 

2.1  Introduction 

The  theory  behind  sampling  plans  which  approximately  minimize  the  maxi¬ 
mum  expected  sample  size  was  presented  by  Weiss  [1]  and  will  be  repeated  here, 
essentially  as  Weiss  presented  it,  with  remarks  relating  to  the  application 
of  the  theory  in  constructing  plans. 

2.2  A  Related  Problem 

The  discussion  will  now  center  on  a  problem  differing  fron  but  related 
to  that  described  by  (1.1).  The  exact  relationship  will  be  de;cribed  later 
in  Section  2.3.  The  problem  to  be  discussed  is  that  of  constricting  a  sam¬ 
pling  plan  which  satisfies  the  conditions 

Pe  (plan  accepts  the  population)  >  1-a 

-1  (2.1) 

(plan  accepts  the  population)  <  6 

and  which  minimizes  the  A.S.N.  curve  at  a  specified  value  0Q(0  i<0o<61)  among 
all  plans  satisfying  (2.1).  It  is  assumed  that  the  observations  taken  are 
independent  Bernoulli  variables.  The  formulation  of  the  above  problem  is 
aimed  at  controling  the  A.S.N.  curve  for  values  of  6  between  and  9j. 

The  following  notation  will  be  employed.  If  T  is  a  plan,  Pq(A|T)  de¬ 
notes  the  probability  that  T  accepts  the  population  when  6  is  the  proportion 
of  defectives.  P0(A|T;Xj , • •  • . Xm)  denotes  the  conditional  probability  that  T 
accepts  the  population,  given  that  T  has  observed  Xj,***,^  and  is  definitely 
going  to  observe  Xm+j  when  0  is  the  proportion  of  defectives.  N  denotes  the 
number  of  observations  that  will  be  taken  before  sampling  is  terminated.  N 
is  a  random  variable  whose  distribution  depends  on  0  and  on  th?  plan  used. 
Ee(N|T)  denotes  the  expected  value  of  N  when  T  is  used  and  0  is  the  proportion 
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of  defectives.  E0 (N |T; Xj, • • •  .X^  denotes  the  conditional  expected  value  of 
N,  given  that  T  has  observed  Xj,***,^  and  is  definitely  going  to  observe 
Xm+1'  w^en  6  *s  t*ie  ProPorti°n  defectives.  Thus  E0(N |T;Xj,  •  •  •  ,\l)  _>m+l. 

If  T  is  a  plan  which  minimizes  E0  (N | T)  among  all  plans  satisfying  (2.1), 
then  there  are  three  positive  constants,  b_j,bQ,bj,  adding  to  unity,  such 
that 

b_1(l-P0  (A | T) ]  ♦  boE0  (N|T)  ♦  bjPe  (A|T)  < 

-1  0  1  (2.2) 

b.Ji-Pg  (A j T * ) )  ♦  boE0  (N|T‘)  ♦  b.(A|T') 

for  each  and  every  plan  T'  (see  [2],  Lemma  4.1).  A  plan  T  satisfying  (2.2) 
it  called  a  Bayes  plan  with  respect  to  the  a  priori  distribution  b_j,b0,bj.  It 
should  be  noted  that  the  Bayes  procedures  are  used  here  only  as  devices  for 
constructing  plans  which  possess  certain  properties.  The  values  of  b_j,bQ,bj 
depend  on  a  and  B,  but  the  exact  relationship  is  unknown.  Therefore,  rather 
than  preassigning  a  and  B,  the  values  b_j.bQ.bj  will  be  preassigned,  and  a 
Bayes  plan  T  will  be  constructed  with  respect  to  b_j,b0,bj.  This  Bayes  plan 
T  then  has  the  following  property.  Jf  T'  is  any  plan  satisfying  the  condi¬ 
tions 

Pe  (A  | T  * )  >  P0  (A  | T) 

-1  -1 

P6i  (A | T ' )  <  P0  ( A | T) 

then  E0  (N | T ' )  >  E0  (N|T).  That  is,  the  values  l-a,B  have  been  replaced  in 
(2.1)  by  P0  ^ (A |T) ,  P0^(A|T),  respectively. 

Define  Q(m,X)  as  equal  to  unity  if  m»0,  and  as  equal  to 

eixi+***4Ni-ei)m-xi . ^ 

e.,Xl*"'*Si-0.1)"-Xl . *■ 

if  m>0.  As  characterized  in  Theorem  4  of  [4],  a  Bayes  plan  T  with  respect 
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to  a  given  a  priori  distribution  b_j,bQ,bj  has  the  following  properties. 
There  is  a  finite  integer  n  and  two  sequences  of  nonnegative  values. 


*01*1  £•••  <*n*b-l/bl>  ro^l  >  •  •  >rn« b.j/bj 
such  that  T  cannot  observe  Xn4^;  and  for  m=0, 1,*  •  •  ,n-l,  T  observes  Xra+1  if 
a^QOn.X)  <  rm>  T  does  not  observe  Xffl+j  and  rejects  the  population  if 
Q(m,X)  >  rm,  T  rrndomizes  in  any  way  between  observing  Xm+1  or  not  observing 
Xfe+1  and  rejecting  the  population  if  Q(m,X)  -rm,  T  does  not  observe  Xm+1  and 
accepts  the  population  if  Q(m,X)  <am,  T  randomizes  in  any  way  between  observ¬ 
ing  Xm+j  or  not  observing  X^+j  and  accepting  the  population  if  Q(m,X)  *3^. 

If  Q(n,X)  =b  j/bj,  T  randomizes  in  any  way  between  accepting  or  rejecting 
the  population.  The  values  n.a^.rg,  all  depend  on  b_i,b0,bi. 

For  the  sake  of  definiteness,  the  only  plans  to  be  considered  will  be 
those  that  definitely  observe  XJn+1  whenever  Q(o,X)  ■  a,,,  or  rm  for  m<n  and 
definitely  accept  the  population  whenever  Q(n,X)  «b_1/b1 .  The  characteriza¬ 
tion  of  T  just  given  is  presented  below. 

n 

If  m«0  then  it  is  understood  that  £  Xj  is  equal  to  zero.  If 

j-1 


m 


*«  '  m 

I  *, 


l  » 

.j-1  V.,) 


m 

-I  Xj 
j-1 


m 


m 

-I  Xj 


<  r_ 


m 


for  o«0, 1, • • •  ,n-l ,  then 

m 


1-01 


_  0i(l-0_i)  - 

108  *»  ‘  108  <ii.id.6D1  *  ■  108  (rreT7)  ‘  rm 


anc  there  is  a  finite  integer  n  and  values  A0*Ai» ’ *  * *An_ j >R0»Ri» *  * ’ *Rn_i ■ 


where 
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,1-01  * 
log  an  -m  logCj^-) 

,  ,0111-0-1*), 
^'eTITW1 


fyn  = 


a-©! 

log  Tm.m  logCjTi-j-) 

.  rei0-8-i), 

l0gt6_1fl-e1)1 


and  A,,,  for  all  m*0,*“,n-l.  For  n 


\  ■  "n 


».l  1-01 

log  gp  -  n  logtpj-j-) 

,  rei(i-e.xJ, 

log[crrnsir1 


(2.3) 


Then  for  m=0, 1, • • • ,n-l,  T  observe.  X  ,  if 

m*  l 

\i  l  *j  i  «m  • 
j-1 

T  does  not  observe  XJn4l  and  accepts  the  population  if 


m 


T  does  not  observe  X,,,^  and  rejects  the  population  if 
m 


T  accepts  the  population  if 
n 


%  • 


If  Ajj,  is  not  an  integer  and  m<n,  can  be  replaced  by  the  smallest  integer 
greater  than  A^,  without  changing  T.  Similarly,  if  Rm  is  not  an  integer  and 
m<n,  Rn  can  be  replaced  by  the  largest  integer  less  than  Rra,  without  chang¬ 
ing  T.  Henceforth,  it  is  assumed  that  the  changes  are  made,  so  that  A^  and 


Rjj  are  integers  for  all  m<n. 


For  a  given  a  priori  distribution  (b.j.bQ.bj),  if  x1#«*  have  been 
observed,  * , X^)  is  defined  as 


bj9jxl*“-,xm  (l-Sj)"-*!-" 


-X, 


m 


_  Xi ♦ • *  ,+X_ 

l  btej  1  ^(l-ej)  1 

i— 1,0,1 


for  j-1,0,1.  The  set  of  quantities  (c.1(fll(Xlf  •••,Xm),  ^(Xj.-.y, 
C1  ra^Xl,”*,XnP)  ca^e^  t*ie  "a  posteriori  distribution  given 
If  m=0  the  quantity  Cj  m(Xj,***,Xm)  *s  defined  as  equal  to  bj .  If  T  is  a 
Bayes  plan  with  respect  to  the  a  priori  distribution  b_j,b0,bi  and  T  has 
observed  Xj,***^^,  then  T  observes  if  and  only  if 


*  co,n,-lCXl-"-'Xm-l)Eeo(NlXl--"-Xm.l) 

*  cl,m-l*Xr"-,Xn-l*Pe1*AlX,'"''Xm-l* 


(2.4) 


*  min^c-ltro-l^Xl,***’Xm-l^cl,m-l^Xl'  *  *  *’Xm-l^  ^ 


for  m*l,**‘,n.  Now  suppose  that  the  values  of  Am  and  ^  are  known  along 
with  the  values  of  P0 _ j (A J Xx , •  •  •  ,Xm) ,  Pei(AlXl’'“*Xm).  and  E0O(NIX1»  *  **»V 
for  all  Xj,***^  with  A^Xj** •  •+Xm<Rp,.  Then  the  following  values  can  be 
computed:  Aml,  l^j,  and  pe  ,  tAlXl‘ ’VP '  P«i(Alxl''"’Xi-l)'  and 

VN|xi’,"’Vi>  for  al1  xr***,xm-i  “ith  ViiV'"‘x.-iiVi-  The 

method  for  computing  the  quantities  mentioned  above  is  presented  in  the  fol¬ 
lowing  two  paragraphs. 

It  is  clear  that  Am  l  ^Ajj,-! ,  for  if  Am  J  <Am-l,  then  if 
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m-J 

.1  xi  *  Vi  • 

i=l 

T  would  observe  and  then  surely  accept  the  population  so  that  observing 
^  would  be  wasteful.  Similarly,  R^  ^  ^R^.  It  should  be  noted  that  for 
most  sampling  plans  of  a  practical  interest  consecutive  values  of  or  Rm 
will  not  be  farther  apart  than  a  value  of  1.  If 
m-1 

i  xi  ■  v1  • 

i=l 


then  for  all  8,  and  m<n, 


petAlxl’ —  ■Xm-1) 

MNIV".Vi> 


i-e+0P0(A|x1,---,xml,i) 

(l-e)m+0Ee(N|X1,---,Xm_1,l). 


(2.5) 


m-1 

I.  xi  ■  R»  • 


i=l 


then  for  all  0,  and  m<n, 


If 


pe(Alxi»,**»xm.i)  ••  (i-e)Pe(A|x1,*'-,xm.i»o) 

E0(Nixr,,,»xm.l)  *  0m+(l-0)E0(N|  Xj,  •  •  •  *Xjn_1»o) 


m-1 

A  -1  <  J  X.  <  R  , 
m  i-1  1  m 


then  for  all  e,  and  m<n, 

VAIV’Vl> 

*  (l-6)Pe(ft|Xj,.. 

E9(M!Xi,...,Xm-i) 

=  (1-0)E0(N|X1,- 


,Xm-l’O^+0P0^AlXl,”',Xm-l’1^ 


xta_1»o)*eEe(N|x1, - •*.xml,i) 


(2.6) 


(2.7) 
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Now  if  inequality  (2.4)  is  satisfied  when 
m-1 

J,  Xi  =  V1  ' 

(using  the  relations  in  equation  (2.5)  to  compute  the  quantities  in  equation 

(2.4))  then  A  .  =  A  -1.  If  (2.4)  is  not  satisfied  when 
m-i  m 

m-1 

l  X.  =  A  -1  , 
i-1  1  m 

then  A  ,  >A  -1,  and  A  ,  is  the  smallest  integer  greater  than  A  -1  for 
m-l  m  m-1  m 

which  (2.4)  is  satisfied  (the  quantities  in  (2.4)  being  computed  by  using 
the  relations  in  the  appropriate  one  of  (2.6)  or  (2.7)).  Similarly,  if  in¬ 
equality  (2.4)  is  satisfied  when 
m-1 

l 


i=  1 


then 


Rffl  If  (2.4)  is  not  satisfied  when 

m-1 

.1  xi  -  R*  • 
i»l 


then  R  .  is  the  largest  integer  less  than  R  for  which  (2.4)  is  satisfied 
m-l  m 

(the  quantities  in  (2.4)  being  computed  by  using  the  relations  in  the  appro¬ 
priate  one  of  (2.5)  or  (2.7)). 

If  the  actual  value  of  n  is  known,  as  for  example  in  Weiss  [10],  ther. 
AjjsRjj  is  known  and  is  given  by  equation  (2.3)  above.  The  quantity  An=Rn  is 
not  necessarily  an  integer  and  cannot  legitimately  be  replaced  by  an  integer 
if  it  is  not  already  one.  Afl_j,  Rn,j»  Pq(a| Xj, •  •  • ,Xn  j)  and  Eq(N|Xj, •• *,X 
can  be  computed  as  described  previously  except  that  on  the  right-hand  sides 
of  (2.5),  (2.6),  and  (2.7) 

MNIXl.***,Xn.i,k)  =  n  for  k=0,l;  and 
Pe(A|Xi.-**.Xn_i,k)  =  1  if  xi+* •  •+Xn_i+k  < An, 

P0(A|X1,***,Xn.1,k)  =  0  if  Xj+- •  •♦Xn.j+k  >  kn,  fork=0,l. 
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Thus,  if  n  is  Known,  the  entire  plan  T  can  be  explicitly  Constructed,  since 

the  values  of  A  ,Aj, • • • ,An,R0,Rj  ,  •  •  •  ,Rn  can  be  found.  Also,  by  repeated  use 

of  the  recursion  formulas  (2.5),  (2.6),  and  (2.7),  Pq(a|T)  and  Eq(N|T)  can 

be  computed  for  any  desired  values  of  6.  The  above  is  accomplished  when  by 

the  recursive  construction  the  point  is  reached  where  the  quantities 

P0(A|no  observations  have  been  taken)  and  Ee(N|no  observations  have  been 

taken)  are  sought.  The  proper  one  of  (2.5),  (2.6),  or  (2.7)  is  then  chosen 

m-1 

by  which  of  the  conditions  the  value  >  Xj  satisfies  (here  this  summation  is 

i-1 

not  precisely  defined  but  is  necessarily  zero  since  no  observations  have  been 
taken) . 

If  the  actual  value  of  n  is  not  known  before  the  construction  of  the 
plan  starts  (as  is  true  for  all  cases  except  the  symmetric  ones  mentioned  in 
Weiss  [10])  but  it  is  known  that  n<.n',  where  n'  is  a  known  finite  integer, 
then  the  Bayes  plan  T  can  still  be  constructed  as  presented  below.  The  con¬ 
struction  technique  mentioned  earlier  is  employed,  proceeding  first  as  though 
n'  were  the  correct  value  for  n.  If  n'  is  actually  greater  than  n,  then  for 
some  value  n*<n',  (2.4)  will  not  be  satisfied  for  any  values  of  xi»*’**xn* 
otherwise,  some  sample  sequences  would  not  terminate  at  n,  which  is  a  con¬ 
tradiction.  This  shows  that  n<^*-l,  so  the  same  construction  is  employed 
again,  proceeding  as  though  n*-l  were  the  correct  value  for  n.  Continuing  in 
this  way,  the  Bayes  plan  T  will  eventually  be  constructed  and  the  value  of  n 
will  be  found  in  the  process.  The  computation  of  n'  will  be  discussed  in 
Section  2.4. 

2.3  The  Relation  Between  The  Two  Problems 

The  Bayes  plan  T  discussed  in  Section  2.2  is  a  generalized  sequential 
probability  ratio  test.  Therefore,  by  Theorem  2  of  Lehmann  [5]  (see  also 
Ghosh  [2]),  if  T  satisfies  (2.1),  it  also  satisfies  (1.1).  This  shows  the 
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relation  between  the  problems  discussed  in  Sections  1  and  2.2,  as  far  as 
the  O.C.  curve  is  concerned. 

It  is  still  necessary,  however,  to  investigate  the  relationship  between 
minimizing  the  A.S.N.  curve  at  a  given  value  90,  and  minimizing  the  maximum 
value  of  the  A.S.N.  curve  for  a  plan  which  satisfies  (2.1)  and  (1.1).  Sup¬ 
pose  T  is  a  Bayes  plan  with  respect  to  the  given  a  priori  distribution 
b_j,bQ,bj.  Let  U  be  any  plan  satisfying  (1.1)  with  a=l~P0  (A|T),  0»P0^(A|T). 
Then  U  satisfies  (2.1)  with  these  values  of  a  and  B.  Define  A(T)  as 

max  E0(N|T)  -  E0  (N|T), 

0  0 

so  that  A(T)  >0.  Then  it  will  be  shown  that 

max  E0  (N|U)  >  max  Ee(N|T)  -  A(T) . 

0  0 

That  is,  T  comes  within  a(T)  of  minimizing  the  maximum  value  of  the  A.S.N. 

curve  among  all  plans  satisfying  (1.1)  with  a=l-P.  (A|T),  0=Pft  (A | T) .  To 

0-1 

show  the  above,  suppose  it  were  not  true.  Then  there  would  be  a  plan  U 
satisfying  the  following  conditions: 

1-P0  x CA|U)  <  1-P0  ^ (a|t) 

p6i(a|u)  <  p0  (a|t) 

Efi  (N|U)  <  max  E.(NjU)  <  max  E_  (N jT)-A(T)  =  Efl  (N|T). 

°o  =  0  0  0  0  o 

But  the  existence  of  a  plan  U  satisfying  these  conditions  would  mean  that 

T  could  not  satisfy  (2.2),  which  is  a  contradiction,  since  T  was  assumed 

to  be  a  Bayes  plan  with  respect  to  b_j,b  ,b^.  This  proves  that 

max  Eq(N|U)  >  max  E  (N |T)  -A (T)  . 
e  0  ~  0  0 

If  A(T)=0,  then  T  actually  minimizes  the  maximum  of  the  A.S.N.  curve. 

In  Weiss  [10]  for  a  special  symmetric  case,  0Q  was  chosen  to  make  A(T)«0. 
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In  general,  however,  there  is  no  direct  means  for  choosing  0Q  exactly;  so  a 
value  of  0O  is  chosen  to  make  A(T)  reasonably  small.  In  Section  2.2  it  was 
pointed  out  that  T  will  continue  sampling  longest  if  Q(m,X)  is  close  to 
b-i/bi  for  all  m.  Define  X(m)  as 

J,  ' 

Then  Q(m,X)  can  be  written  as 


el(l-0_l) 

mX(m) 

1-0! 

e.jd-ep 

i-0-! 

If  Q(m,X)  is  close  to  '>_j/bi,  then 

vlX(m) 


e1(l-9.1) 

O.i(l-Qi) 


is  close  to 


1-M 

1 

»— » 

_ k 

l1-0!  J 

L*>i  J 

il/m 


which  for  large  m  is  close  to  (1-8_j)/(1-6j)  if  b.^/bj  is  not  too  far  from 
unity.  However,  for  large  m,  X(m)  is  with  high  probability  close  to  e,  the 
true  proportion  of  defectives  in  the  population.  Then  it  would  be  expected 
that  the  A.S.N.  curve  is  high  at  a  value  of  e  such  that 


'0lU-0_l) 

0_i(i-0iV 


e 


is  close  to 


1-0 

1~ 


or  at  a  value  of  6  close  to 


log[ 


1-0 

He 


-1 

1 


] 


©id-Q-i), 

0.l(l-0!)J 


(2.8) 


If  b.j/bj  is  not  close  to  unity  and  if  there  is  approximate  information 
available  as  to  the  largest  number  of  observations  the  test  could  actually 
take  (perhaps  by  interpolation  among  the  Bayes  plans  presented  later),  then 
the  maximizing  0  will  be  approximately,  for  large  n+, 
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-1  1  D-l 

“^rep  *  ?  108 ^T1 
e,(i-e  .) 


(2.9) 


where  represents  the  approximate  information.  If  the  A.S.N.  curve  is 
minimized  at  the  value  of  0  given  by  (2.8)  or  (2.9),  it  may  be  expected  that 
the  maximum  value  of  the  A.S.N.  curve  will  be  approximately  minimized. 
Therefore  the  value  of  either  (2.8)  or  (2.9)  is  used  in  constructing  the 
Bayes  plan  T,  with  the  expectation  that  this  choice  will  make  A(T)  small. 


jer  Bound  on  the  Largest  Possible  Sample  Size 


The  computation  of  the  quantity  n',  an  upper  bound  for  n,  will  be  dis 


cussed  here.  The  quantity  n'  is  a  function  of  b_j.b0,bj,  but  for  typograph¬ 
ical  simplicity,  the  notation  will  not  exhibit  this  dependence. 


It  would  be  ideal  if  the  exact  value  of  n  were  available  before  start¬ 


ing  the  construction  of  the  Bayes  plan  T.  However,  at  present  the  exact 
value  can  be  found  only  in  special  symmetric  cases  (as  in  [10]);  so  that  n' 
must  be  used,  relying  on  the  actual  construction  of  T  to  find  the  value  of 
n.  The  cruder  (greater)  the  upper  bound  n',  the  greater  the  amount  of  com¬ 
puter  time  wasted  in  working  back  to  n.  Below,  a  value  of  n'  is  developed 
which  is  used  in  the  construction  of  T.  Wetherill  [11]  and  Ray  [6]  have 
studied  a  different  upper  bound  on  n  which  has  not  been  compared  to  the  bound 
presented  here. 

From  the  discussion  immediately  following  Lemma  4.1  and  in  Lemmas  4.3 
to  4.6  in  [4],  it  is  found  that  there  is  a  continuous  nonnegative  function 
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T  does  not  observe  X  ,  if 

m*i 

c,  m(X  *..,X  ) 

c  (X,,...,X  )  >  K[— - JL]  ; 

°'m  1  "  c4.«lXl'-"'V 

T  can  randomize  in  any  way  between  observing  or  not  observing  Xm+j  if 

c  (X  ,...,X  ) 

c  (X,,...,X  )  =  K[  ; 

O'”  1  ”  cl,m(Xl'  —  'V 

and  if  T  stops  sampling  immediately  after  observing  X^,  T  accepts  the  popul? 
tion  if 


C4,n,<Xl'-"'Xm' 


<  1  , 


T  rejects  the  population  if 


ci,»cV"-V 


>  i . 


T  randomizes  in  any  way  between  accepting  or  rejecting  if 


Ci,»(xi-->V 


L.  I  A,  ,  •  •  •  ,  A  1 

4,mv  1’  *  nr 


=  1  . 


It  should  be  noted  that  the  function  K(u)  does  not  depend  on  b_j,b0>bj. 
Furthermore,  by  the  convexity  of  the  acceptance  and  rejection  regions  de¬ 
scribed  by  Kiefer  and  Weiss  in  [4],  K(u)  has  the  following  properties: 

max  K(u)  =  K(l)  <  1  ,  K(0)  =  K(»)  =  0, 

u>p 

and  if  u^.Uj  are  any  two  values  satisfying  one  of  the  conditions  (teu^u^l 
or  l<ui<U2,  then  for  all  d  in  the  closed  interval  [0,1] 


1-K(u,)  l-K(u-) 

dUjI-j— — )  .  (l-dJUjf-j— ) 

X  I - itktsts - rTTuT-1  <  dKfup.Cl-dJKtup 


(2.10) 
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In  particular,  setting  u1»0,u2=0  in  (2.10)  and  denoting  K(l)  by  k,  it  is 
found  that 

„  r  H(l-'i)  (1-k)  ,  .  „ 

K  =  (1‘d)k 

for  all  d  in  [0,1].  Setting  Uj=l,u2*“  in  (2.10),  it  is  found  that 


(2.11) 


*  Ad(l-k)t(l-d),  .  Jt. 
K  [  ~  ijdcr-vi — ]  - dk 


(2.12) 


for  all  d  in  [0, 1] . 

If  the  function  K(u)  were  known,  the  entire  test  T  could  be  constructed 
explicitly,  and  the  exact  value  of  n  would  be  known.  However,  K(u)  is  not 
known.  Suppose  that  a  known  function  L(u)  satisfies  L(u)  ^K(u)  for  ill  u>0, 
and  that  n'  is  defined  as  the  largest  possible  sample  size  of  a  plan  con¬ 
structed  by  acting  as  though  K(u)  were  equal  to  L(u) .  Then  it  is  clear  that 
n'>n.  Furthermore,  the  closer  L(u)  is  to  K(u),  the  smaller  n'  will  be.  The 
next  step  is  to  construct  a  function  L(u)  which  is  known  to  satisfy  L(u)  ^K(u) . 

It  can  be  said  that  a  nonnegative  value  u  satisfies  the  condition  S  if 
no  Bayes  plan  with  respect  to  the  a  priori  distribution 


^  ■  ««>  • 

can  ever  observe  X2.  From  the  characterization  above,  there  is  a  Bayes  plan 
TQ  with  respect  to  the  a  priori  distribution  which  does  not  observe  Xj,  and 
there  is  a  Bayes  plan  T  which  does  observe  X^.  Since  E0  (n|tq)=0  and 
P0j(A|To)»Pq  (A|Tq),  Tq  makes  the  left-hand  side  of  (2.2)  equal  to 

,1-K(u)  rl-K(u)11 

B,n{T*u  - 

Since  Tj  observes  X^  but  cannot  observe  X2,  it  is  clear  that  Tj  accepts  the 
population  if  X^=0  and  rejects  the  population  if  Xj*l,  otherwise  it  would 
not  be  necessary  for  T^  to  observe  X^.  Then  E0  (N|T,)=1,  P0  ^  f  A I T*  1 »-  ■*  *.i. 
Pe  (AlT1)-l-e1.  so  that  Tj  mnkca  o.v  ».<*  w—J  cide  of  (2.2)  equal  to 
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Since  both  TQ  and  Tj  minimize  the  left-hand  side  of  (2.2),  then 

rl-K(u),  w  .  rl-K(u)i  •  / l-K(u)  ,1-K(u) , , 

Using  this  equality,  it  is  found  that  if  u  satisfies  the  condition  S,  then 


l-K(u) 


K(u)  = 


uere.i 

1+U+U0J-0J 


if  u  <  1 


(2.13) 


1-U+U0.-0 

-  ;«ue,:e'  if  “i1  • 


Next  it  will  be  shown  that  u=l  satisfies  the  condition  S,  so  that 


k  =  K(l)  = 


V6-! 

2+0  -e  ’ 


Using  the  a  priori  distribution  %(l-k) ,k,%(l-k) ,  it  is  found  that  if  Xj*0, 


k^-v  ci,i(xi)  ^i 

co,ilV  =  k(i-e0)**s(i-k> (2-e1-e_1)  '  c4(1(Jf  }  “  i-e_1 


SU-k^-e^) 


V°-i- 


in  (2.11),  it  is  found  that 


1  -  6  j  kd-ep 

K[rr[]  -  i-61^(i-k)(2-er0_1] 


and  it  is  easily  verified  that 


kd-ep 


k(l-0o) 


( e  i  -  e  _  x ; 


Therefore  when  X^=0, 


C1  1^1^ 
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T1 


so  X2  cannot  be  observed.  If  X^»I, 

ke„ 

C  ,(X,J 


C1,1(V 


o,iui  3  ke0^d-k) (el+e_1)  *  '  e 


-1 


if 


d  «  1 


*s(i-k)  fe1-e_1) 


e.i+1sCi-k)  Ce1-e_1) 


in  (2.12),  it  is  found  that 


6j  ke 

K  'r^  i  e^vqrr-TTTe'p.jj  * 


and  it  is  easily  verified  that 


k8 


-1 


ke. 


0_1+is(l-k)(e1-e_1)  <  keo*Jj(i-kj  (ej+e.j)  * 


Therefore  when  Xj*l, 

co,l(Xl)  >  X  • 


cl,ltxl) 


so  X2  cannot  be  observed.  This  proves  that  u=l  satisfies  the  condition  S. 


Now,  a  function  Lj(u)  is  defined  as  follows: 


Ljfu)  = 


l^(l-k)(~) 


if  u  <  1  , 


Ll(u)  *  l+4(l-k) (u-1) 


if  u  >  1  . 


Using  (2.11)  and  (2.12),  it  can  be  verified  that  Lj(u)  ,>K(u)  for  all  uj>0. 


Also  Lj(u)  is  a  known  function  of  u,  since  k*K(l)  has  been  computed  explic¬ 


itly  above.  Thus  Lj(u)  could  be  used  as  the  function  L(u)  described  above. 
However,  more  can  be  done  by  finding  values  of  u  other  than  unity  which 
satisfy  the  condition  S. 


J. 
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and  it  is  evident  that  l(u)  ^KCu)  for  all  u^o.  Now  the  explicit  values 
for  v  and  w  can  be  computed. 

First  it  is  noted  that  K(u)  satisfies  the  following  inequalities: 


K(u) 

.  uVe_i 

“  l+u*u6^-6  j 

if  u<l 

(2.15) 

K(u) 

1-U4U0J-6  j 

if  u  >  1  . 

_ac 

Note  that  the  right-hand  sides  of  these  inequalities  are  identical  with  the 
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l 


right-hand  sides  of  the  equalities  (2.13)  which  hold  only  for  values  of  u 
satisfying  the  condition  S,  whereas  (2.15)  holds  for  all  u.  The  proof  of 
(2.15)  proceeds  similarly  to  that  of  (2.13),  noting  that  TQ  is  a  Bayes  plan 
with  respect  to  the  a  priori  distribution 


l-K(u) 

1-u 


K(u) 


u  [ 


l-K(u), 
1+u  J 


But  Tj  may  not  hold  for  general  u,  and  therefore  the  left-hand  side  of  (2.2) 
using  Tq  is  less  than  or  equal  to  the  left-hand  side  of  (2.2)  using  Tj .  The 
inequalities  (2.15)  then  follow  directly. 

Define 


0 


max  { 


1 


e.+e  ,-e  e,-e  ,e, 

1  -1  Ol  -11 


0‘1  + 


9 


1 


Then  o<Vj<l,  and  a  straightforward  but  somewhat  tedious  calculation 
shows  that  if  u  is  the  open  interval  (v^l)  and  if  a  Bayes  plan  is  constructed 
with  respect  to  the  a  priori  distribution 


l-K(u) 


,  K(u)  ,  u  [ 


l-K(u) i 
1+u  J 


» 


then 


cm'V 


co,l(Xl)  *  L1  (c  jjtXj)1 


For  Xj=0  or  1,  (the  first  inequality  in  (2.15)  is  used  in  the  calculation). 
Since  Lj(u)  ^Kfa)  for  all  u,  then  any  u  in  (Vj,l)  satisfies  the  condition  S. 
t^fine 


mi-  { 


(l-0o)(l-0.j) 

a-01)(ite1-eo-e_1)  * 


1*e-r(5T':)(ei 


■e-P 


i-e, 


Then  Wj>1,  and  calculation  stw,  that  if  u  is  in  the  open  interval  (l,w^) 
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and  if  a  Bayes  plan  is  constructed  with  respect  to  the  a  priori  distribution 


l-K(u) 

“T+u 


K(u) 


then 


ci  i(xi) 

C°»  1  (Xl^  ”  Ll  t 


For  XjcO  or  1,  the  second  inequality  in  (2.15)  is  used  in  the  calculation. 
Since  L^fu)  >^K(u)  for  all  u,  then  any  u  in  (l,w  )  satisfies  the  condition  S. 
Now  define  L2(u)  as  L(u)  was  defined  in  (2.14),  using  v^.w^  in  place  of  v,w. 
Clearly,  L2(u)^K(u)  for  all  u.  Noting  that  v^,w1  is  found  by  using  the  fact 
that  Lj(u)  >^K(u)  for  all  u,  and  noting  that  Lj(u)  ^L2(u)  ^K^)  for  all  u,  it 
would  be  possible  to  find  values  v2,w2,  such  that  all  u  in  (v2,w2)  satisfy 
the  condition  S,  where  v2<vi»w2>wi’  The  quantities  v2,w2  would  be  computed 
using  L2(u)  just  as  Lj(u)  was  used  in  computing  Vj,Wj.  Then,  using  v2,w2  in 
(2.14),  there  would  be  a  function  L2(u)  with  L2(u)  >_L2(u)  ^K(u) .  The  above 
process  could  be  continued.  Since  the  establishment  of  an  upper  bound  n'  may 
be  of  only  relative  importance  in  the  wholesale  construction  of  the  Bayes 
plans,  the  known  function  Lj(u)  defined  by  using  v,w  in  (2.14)  was  used  here 
as  the  final  function  L(u). 

For  a  given  a  priori  distribution  b.j,b0,bj  define  z.m)  as  the  largest 
integer  for  which  Q(m,X)  ^b.j/bj  when  Xj+* •  •+Xn  =  z(m)  for  each  positive  in¬ 
teger  m.  Then  from  Section  2.2,  it  is  known  that  if  a  Bayes  plan  cannot  ob¬ 
serve  Xm+^  when 

m 

l  X.  =  z(m) 
i=l 


or  z(m)+l,  the  plan  cannot  observe  X  j  under  any  circumstances.  Define  c(m) 
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as  the  value  of  c _ (X,,*««,X  )  when 

o,m  l  in 

m 

l  x.  =  z(m)  , 

X  _  1  A 


and  define  c'(m)  as  the  value  of  cq  w^en 

m 

[  X.  *  z(m)*l  . 
i*l  1 


Define  f(m)  as  the  value  of 


‘l.m'V-V  . 


and  f'(m)  as  the  value  of 


when 


m 

I  xi 

i-i 


m 

l  h 


~  z(m)  , 


z(m)+l  . 


Define  n'  as  the  smallest  integer  satisfying  both  of  the  following  inequali¬ 
ties: 


c(n')  >  Lj(f (n1)) 
c'(n')  >  L^f'Cn'))  . 

If  Lj(u)  were  equal  to  K(u),  the  Bayes  plan  could  not  observe  Xn,+j.  Since 

it  is  known  that  Lj(u)  ^K(u)  for  all  u,  it  is  also  known  that  the  Bayes  plan 

can  never  observe  X  ,  . .  Therefore  n'  is  the  upper  bound  for  n  described  in 

n'+i  rr 

Section  2.2. 

Since  n'  is  a  discrete  variable  and  bQ  is  a  continuous  variable,  in  the 
actual  computation  it  is  easier  to  reverse  the  above  process  and  search  for 
bQ,  rather  than  for  n'.  The  reverse  process  is  completed  by  fixing  a  posi¬ 
tive  integer  m  and  a  value  for  the  ratio  b.j/bjj  then  z(m),  f(m),  and  f'(ni) 
can  be  computed  without  knowing  the  value  of  bQ.  Define  b^(m,  b.j/bj)  as 
the  smallest  value  of  bQ  that  makes  both  of  the  following  inequalities  hold: 
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c(m)  >,  L^ffm)) 
c'(m)  >  Lj (f ' (m) )  . 

Then  it  is  clear  that  there  is  a  Bayes  plan  with  respect  to  the  a  priori  dis¬ 
tribution  b  j.b^Cm.b  j/bj)  .bj  which  never  observes  where  in  this  a 

priori  distribution  b  j/bj  is  equal  to  the  fixed  given  ratio.  Then  the  en¬ 
tire  Bayes  plan  can  be  constructed  by  working  back  from  m  as  the  upper  bound 
on  the  number  of  observations. 


I 


' 

1 

I 


3.  Sequential  Sampling  Plans  Which  Approximately  Minimize  the  Maximum 
Expected  Sample  Size 

The  format  of  each  plan  (the  order  of  plans  is  listed  in  the  table  of 
contents  by  values  of  0_p6^,oi,B)  is  similar  to  the  following  example: 


rH 

o 

o 

II 

H 

1 

<x> 

9  j  =  .011  a  =  | 

8  =  .01 

Comparison  of  A.S.N. 

Values 

Test 

6  =  .001 

Maximizing  e 

Minimax 

611.5 

735.0 

Wald 

568.9 

885.0 

Fixed  Sample  Size 

965.0 

Plan  1 


e  =  .oil 

321.0 

28S.8 


n 

An 

Rn 

1 

n 

Rn 

1 
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The  Minimax  test  refers  to  the  one  introduced  by  Weiss  [1].  A.S.N.  values 
for  the  Wald  test  [9]  were  computed  using  the  Wald  approximation  formulas, 
and  the  number  of  observations  required  by  the  fixed  sample  size  test  was 
computed  using  a  normal  approximation.  The  a, 6  values  are  approximately  those 
for  the  given  plan.  The  exact  a, 6  values  can  be  found  in  [3].  A.S.N.  values 
for  the  Minimax  and  Wald  tests  are  compared  at  0 _ j and  the  "Maximizing  9." 
The  value  of  6  associated  with  the  maximum  A.S.N.  is  usually  different  for 
the  two  sequential  tests. 

The  acceptance  and  rejection  numbers  at  stage  n,  are  labeled  An,Rn  re¬ 
spectively.  The  batch  of  product  being  sampled  is  accepted  if  the  cumulative 
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number  of  defectives  after  the  n—  observation  is  no  greater  than  and  is 
rejected  if  the  cumulative  number  of  defectives  is  no  less  than  R^  Other¬ 
wise,  another  observation  is  taken,  and  the  comparison  is  repeated.  If  the 
Ah  column  is  dashed  for  the  n—  observation,  the  batch  cannot  yet  be  accepted, 
and  if  the  Rn  column  is  dashed  for  the  n—  observation,  the  batch  cannot  yet 
be  rejected.  The  values  of  A^  and  Rn  are  only  listed  when  one  or  both  change 

value.  For  example,  in  Plan  1  presented  above  no  decision  to  accept  or  reject 

Tci 

can  be  made  until  three  observations  have  been  taken,  for  the  3 —  through  the 
s  t 

241 —  observations  the  batch  cannot  be  accepted  but  can  be  rejected  if  a  total 
of  3  or  more  defectives  is  observed,  etc.  The  last  entries  in  the  Ah  and  Rh 
columns  differ  by  one  unit  and  force  sampling  to  terminate  if  it  reaches  this 
point. 
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13  ABSTRACT 


The  author,  in  his  masters  thesis,  constructed  a  catalog  of  442 
Bernoulli  sampling  plans  which  approximately  minimize  the  maximum  expected 
sample  size  among  all  plans  which  guarantee  certain  O.C.  probability  require¬ 
ments.  Fifty-two  of  these  plans  (which  would  appear  to  be  of  greatest  prac¬ 
tical  interest)  are  presented  in  this  report.  A.S.N.  curve  comparisons  are 
made  with  plans  based  on  the  Wald  sequential  probability  ratio  test  and  the 
fixed  sample  size  test  which  guarantee  the  same  O.C.  probability  requirements, 
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